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Abstract 



6 '• 

p! Context. ChandrafHETG observations of the recurrent nova RS Ophiuchi at day 13.9 of its 2 006 outburst reveal a spectrum coveri ng a large 
^ range in plasma temperature and characterized by asymmetric and blue-shifted emission lines dNelson et al.ll2008l : lDrake et alj2008h . 

Aims. We investigate the origin of asymmetries and broadening of the emission lines observed with ChandrafHETG. We explore possible 
diagnostics of the early blast wave and of the circumstellar medium (CSM) in which the explosion occurred. 

Methods. We perform 3-D hydrodynamic simulations of the blast wave from the 2006 outburst, propagating through the inhomogeneous 
CSM. The model takes into account the thermal conduction (including the effects of heat flux saturation) and the radiative cooling. From the 
simulations, we synthesize the X-ray emission and derive the spectra as they would be observed with ChandrafHETG. 
• Results. The simulated nova remnant is highly aspherical and the blast wave is efficiently collimated by the inhomogeneous CSM. Our model 
reproduces the observed X-ray emission in a natural way if the CSM in which the outburst occurred is characterized by an equatorial density 
enhancement. Our "best-fit" model predicts that most of the early X-ray emission originates from a small region propagating in the direction 
t- —i • perpendicular to the line-of-sight and localized just behind the interaction front between the blast wave and the equatorial density enhancement. 
"— "i , The model predicts asymmetric and blue-shifted line profiles remarkably similar to those observed. These asymmetries are due to substantial 
■ X-ray absorption of red-shifted emission by ejecta material. 
Conclusions. The comparison of high quality data of Chandra/HETG with detailed hydrodynamic modeling has allowed us to unveil, for the 
first time, the details of the structure emitting in the X-ray band in early phases of the outburst evolution, contributing to a better understanding 
? "~j ■ of the physics of interactions between nova blasts and CSM in recurrent novae. This may have implications for whether or not RS Ophiuchi is 
K\ a Type la SN progenitor system. 
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1. Introduction 

RS Ophiuchi (RS Oph) is a symbiotic recurrent nova that 
went into its latest outburst on 2006 February 12.83 UT 



dNarumi et al 



1898, 1933, 



2006); previous outbursts we re recorded in 
1958, 1967, and 1985^ (see iRosinol Il987t 



Rosino & Iiiimal 119871 and references therein). RS Oph is 
thought to be a binary system, comprising a red giant star that 
does not fill its Roche lobe, and a white dwarf of mass near the 
Chandrasekhar limit ( Dobrzvcka & Kenv on 1994: IShore et al 



19961: iFekel et al.ll200(i; the system has a period of 455.72 
0.83 days dFekel et al.ll2000l) . Current theory tells us that, in 
this class of binaries, the outbursts occur on the white dwarf 
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1 T wo possibl e outbursts were also observ ed in 1907 dSchaefeJ 
l2004h and 1945 foppenheimer & Matteil 1993b . 



due to thermonuclear runaway of hydrogen-rich material ttans- 
ferred from the companion red giant onto the surface of the 
white dwarf (e.g. lKatolll99lt lAnupama & Mikolaie wskall 19991: 
Starrfield et al.ll2000l) . 



The strong interest in studying RS Oph originates mainly 
from the belief that recurrent novae are progenitors of Type la 
supernovae (SNe la). Two possible scenarios have been pro- 
posed: i) matter from the red giant is accreted by the white 
dwarf, causing the latter to increase in mass until it reaches the 
Chandrasekhar limit and explodes as an SN la; ii) two white 
dwarf stars merge, with the combined mass momentarily ex- 
ceeding the Chandrasekhar limit, causing an SN la explosion. 
Which of these two scenarios is the more plausible is still un- 
der debate, so studying the nova outbursts from RS Oph can 
provide important insight into SNe la. 

During the 2006 outburst an intensive international observ- 
ing campaign was organized, incorporating observations rang- 
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ing from radio to X-ray wavelengths, and monitoring the out- 
burst since the early phases of its evolution. The results in- 
cluded the detection of X-ray emission from hot gas, evolving 
from a dominan t plasma temperature o f ~ 10 KeV few days af- 
ter the eruption Jsokoloski et al.ll2006l) to » 4 KeV 10 days af - 
ter optical maximum (IBode et al . 2006; Sokoloski et al . 2006). 
It has been suggested that bright X-ray emission arises from 
the shock-heated extended outer atmosphere of the red giant 
(e.g.lBode&Kahnl l 19851: IBode et alJl20o5 iNelson et alJl2008 ; 



Drake et alj|2008l) . 

Chandra/HETG observations at day 13.9 revealed a rich 
spectrum of emission lines indicative of emitting plasma with 
temperatures ranging be t ween 3 and 60 M K (INelson et al 



20081: brake et al.1 120081) . brake et all d2008l) noted that the 
lines are too strongly peaked to be explained by a spherically- 
symmetric shock, suggesting a collimation mechanism of the 
X-ray emitting plasma in the direction perpendicular to the 
line-of-sight. The lines also appear asymmetric and slightly 
blue-shifted, while the red wings of the line pro fi les be come 



weaker with increasing wavelength, brake et al.l (120081) sug 



gested that the asymmetric nature of the circumstellar medium 
(CSM) in which the explosion occurred can be responsible for 
both the broad range in plasma temperature and the shock col- 
limation observed. 

Here we investigate the origin of the line asymmetries, 
broadening and blue-shifts observed with Chandra/HETG . We 
aim at exploring possible diagnostics of the early blast wave 
and of the inhomogeneous CSM in which the explosion oc- 
curred. To this end, we model the expansion of the blast wave 
from the 2006 outburst of the recurrent nova RS Oph through 
the extended outer atmosphere of the companion red giant, us- 
ing detailed 3-D hydrodynamic simulations. From the simula- 
tions we synthesize the X-ray emission and derive the spectra 
as they would be observed with Chandra/HETG. 

In Sect.|2]we describe the hydrodynamic model, the numer- 
ical setup, and the synthesis of X-ray emission; in Sect. [3] we 
discuss the results; and finally in Sect. [4] we draw our conclu- 



2. Model 

2.1. Hydrodynamic modeling 

The blast wave is modeled by numerically solving the time- 
dependent fluid equations of mass, momentum, and energy 
conservation in a 3-D Cartesian coordinate system (x, y, z), tak- 
ing into account the radiative losses from an optically thin 
plasma and thermal conduction (including the effects of heat 
flux saturation): 

-| + V.pu 



. 



— + V • puu + VP = , 

dt H 

— + V • (pE + P)u = -V • q - n e n H A(T) 
dt 



(1) 



is the total gas energy (internal energy, e, and kinetic energy), t 
is the time, p - pm H n^ is the mass density, p = 1 .3 is the mean 
atomic mass (assuming cosmic abundances), m H is the mass of 
the hydrogen atom, «h is the hydrogen number density, n e is the 
electron number density, u is the gas velocity, T is the temper- 
ature, q is the conductive flux, and A(T) r epresents the radia- 
tive losses per unit emis s ion measure (e.g. Raymond & S mith 
11977k IMewe et al.lll985b IKaastra & MewdboOOh . We use the 
ideal gas law, P - (y - l)pe, where y = 5/3 is the adiabatic 
index. 

To allow for a smooth transition between the classical and 



saturated conduction regime, we followed iDalton & Balbus 
d 19931) and defin ed the conductive flux as (see also 



Orlando et alj|2005h 



1 1 
2 = 1 — + — 

\*7spi <?sat 



(2) 



Here q sp i represents the classical conductive flux (Spitzer 1962) 



-6 T e20\- 



3/2 



(kBT) 5/2 k b 
me /2 e 4 Zln(A) 



vr 



(3) 



where is the Boltzmann constant, m e is the electron mass, 
e is the electron charge, Z is the average atomic number, 
ln(A) is the Coulomb logarithm, 6j and e parameters de- 
pend on the chemical composition, and in a proton-electron 
plas ma: 6t - 0-225 and e = 0.419. The saturated flux, g sat , 
is (ICowie & McKedl 19771) 



9sat 



-sign(Vr) 5<ppc s , 



(4) 



where c s is the isothermal sound speed, and <p is a number on 
the order of unity. We set <f> = 0.3 according to the va lues sug- 
geste d for a fully ionized cosmic gas: 0.24 < <fi < 0.3 5 dGiuliani 



1984; Borkowski et al. 1989, Fadevev et al. 2002, and refer- 



ences therein). 

In order to trace the motion of the material ejected by the 
violent eruption, we consider a passive tracer associated with 
the ejecta. To this end, we add the equation 



dt 



+ V • C cl u = 



(5) 



Here 



E = e+ -Ini- 



tio the standard set of hydrodynamic equations. C c i is the mass 
fraction of the ejecta inside the computational cell. The ejecta 
material is initialized with C c i = 1, while C c i = in the ambient 
medium. During the shock evolution, the ejecta and the ambi- 
ent medium mix together, leading to regions with < C c \ < 1 . 
At any time t the density of ejecta material in a fluid cell is 
given by p d =pC d . 

The calculations described in this paper were performed 
using flash, an ada ptive mesh refinement multiphysics code 
(IFrvxell et al. I l2000h . The hydrodynamic equations are solved 
using th e flash implementation of the piecewice-parabolic 
method (IColella & Woodward! 1 19841) . The code was designed 
to make efficient use of massively parallel computers using the 
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message-passing interface (MPI) for interprocessor communi- 
cations. The code has been extended with additional computa- 
tional mo dules to handle rad iative losses and thermal conduc- 
tion (see lOrlando et alj |2005 for the details of the implementa- 
tion). 

As an initial condition, we assume a spherical Sedov-type 
blast wave with radius ru) = 1/3 AU and with total energy 
£bo, originating from the thermonuclear explosion on the white 
dwarf. The initial total energy of the explosion is partitioned 
so that 1/4 of the energy is contained in thermal energy, and 
the other 3/4 in kinetic energy. The initial total mass of the 
ejecta is M e j. The blast propagates through the extended outer 
atmosphere (the wind) of the companion red giant and is off-set 
from the origin of the wind d ensity distribution by 1 .5 AU (i.e. 
the sys tem orbital separation ; lDobrzvcka& Kenvonll994 .We 
follow lO'Brien et alj ( 12006b and assume the gas density in the 
red giant wind proportional to r~ 2 (where r is the radial distance 
from the giant) and its temperature 2 x 10 4 K. Note that the 
temperature values of the red giant wind (in symbiotic stars), 
r w , ca n range between few 10 3 K up to 10 5 K closer to the red 



giant (Friedjung et al. 1983). To evaluate the effects of T w on 



the results of our simulations, we compared simulations with 
T„ = 4000 K with simulations with T w = 20000 K and found 
that the results do not change. 

In addition to the r~ 2 density distribution, we include 
also an equatorial density enhancement (hereafter EDE) in 
the red giant wind, as suggested by VLBA radio synchrotron 
observation s ( O'Brien et al.1 120061) and by HST observations 
dBode etalJl2007h of the 2006 blast wave, and as predicted 
by detailed hydrodynam ic modeling (Mas trodemos & Morris 
1999tlWalder et al.ll2008l) . The mass density distribution of the 
unperturbed CSM is given by: 



Pw 



•Peq ' 



-(/i/L) 2 



(6) 



where p w = [rniun^, is the mass density at a distance of 1 AU 
from the red giant, r au is the radial distance from the giant in 
AU, peq = ^'«H«eq is the density enhancement at the equato- 
rial plane, h is the height above the equatorial plane and L is a 
characteristic length scale. 

In our 3-D simulations, the symmetry of the problem al- 
lows us to solve the hydrodynamic equations in one quadrant 
of the spatial domain (see Fig . |T} . The coordinate system is ori- 
ented in such a way that both the white dwarf and the red giant 
lie on the z axis. The red giant is at the origin of the coordinate 
system, (x, y, z) = (0, 0, 0), and the computational domain ex- 
tends 60 AU in the z direction, and 30 AU in both the x and y 
directions; the white dwarf is arbitrarily located to the north on 
the z axis (x = y — 0) at z = 1.5 AU (i.e. the system orbital 
separation; Dobrzycka & Kenyonlll994 ). 

At the coarsest resolution, th e adaptive mesh algo rithm 



used in the flash code (paramesh; MacNeice et al. 2000) 



uni- 
formly covers the 3-D computational domain with a mesh of 
4x4x8 blocks, each with 8 3 cells. We allow for 8 nested 
levels of refinement during the first 3 days of evolution and 
for 6 levels in the rest of the simulation, with resolution in- 
creasing twice at each refine ment level. The refinement cri- 
terion adopted (ILohnerll 1987b follows the changes in density, 



white dwarf 



red giant 




Figure 1. Sketch of the initial geometry of the RS Oph sys- 
tem. The figure shows the quadrant of the spatial domain that 
is modeled numerically. The red giant is at the origin of the 
coordinate system, (x,y,z) - (0,0,0), and the computational 
domain extends 60 AU in the z direction, and 30 AU in both 
the x and y directions; the white dwarf is located on the z axis 
(x — y = 0) at z — 1.5 AU. The (y,z) plane is the equato- 
rial plane and the green structure marks the equatorial density 
enhancement. 



temperature, and tracer of ejecta. This grid configuration yields 
an effective resolution of « 0.03 AU at the finest level dur- 
ing the first 3 days of evolution (when the radius of the nova 
remnant was < 5 AU) and as 0.1 AU at later times, corre- 
sponding to an equivalent uniform mesh of 1024 x 1024 x 2048 
and 256 x 256 x 512 grid points, respectively. The higher spa- 
tial resolution during the first 3 days of evolution is required 
to describe properly the development of initially small struc- 
tures dominated by radiative cooling. We use reflecting bound- 
ary conditions at x = x m m and y - y m ; n (consistent with the 
adopted symmetry) and zero-gradient conditions at the other 
boundaries. 

We follow the expansion of the blast wave through the red 
giant wind for « 15 days, considering two sets of simula- 
tions: with or without the EDE. In both cases, we explore a 
parameter space defined by the initial mass of ejecta M e j in the 
range 10~ 7 - 10~ 6 M , the initial energy of the explosion Eq in 
the range 10 43 - 10 44 erg (namely the published e stimates for 



the ejected mass and energy of the explosion, e.g . lYaron et al 
20051 ISokoloski et ai]|2006t b'Brien et ai]|2006h . and the par- 
ticle number density of the red giant wind at 1 AU n w in the 



range2xl0 7 -2xlO I0 cm 3 (e.g. lDobrzvcka & Ken vonl 19941) . 



In the set of simulations with the EDE, we set n en = 10 cm 



(consistent with the values found bv lWalder et al.l2008l) and we 
explore the case of L — 1 and 2 AU (see Eq.|6]l. Table[T]summa- 
rizes the physical parameters characterizing the simulations. 
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Table 1. Initial parameters of the 3-D simulations 



the value determined from HI 21 cm measurements, e.g. by 
Hjellming et al . 1986, an d consistent with a distance to RS Oph 



run 


M". 


F 




A 


u 




[M G ] 


[erg] 


[cm' 3 ] 


[cm- 3 ] 


[AU] 


ND-E43-N7 


10-' 


10 43 


2 x 10 7 







ND-E44-N7 


lfr 6 


10 44 


2 x 10 7 







ND-E43-N8 


io- 7 


10 43 


2x 10 8 







ND-E44-N8 


lfr 6 


10 44 


2x 10 8 







ND-E43-N10 


lfr 7 


10 43 


2x 10 10 







ND-E44-N10 


lfr 6 


10 44 


2x 10 10 







YD-E43-N7-L2 


lfr 7 


10 43 


2x 10 7 


10 s 


2 


YD-E44-N7-L2 


lfr 6 


10 44 


2 x 10 7 


10 s 


2 


YD-E43-N8-L2 


lfr 7 


10 43 


2 x 10 8 


If) 8 


2 


YD-E44-N8-L2 


lfr 6 


10 44 


2 x 10 8 


10 s 


2 


YD-E44-N7-L1 


lfr 6 


10 44 


2x 10 7 


10 s 


1 



" Initial mass of ejecta; 

b Initial energy of the explosion; 

c Particle number density of the red giant wind at 1 AU; 

d Particle number density enhancement at the equatorial plane; 

e Characteristic length scale in Eq.|6) 



2.2. Synthesis of the X-ray emission 

From the model results, we synthesize the X-ray emission orig- 
inating from the blast w ave, applying a meth odology analogous 



to the one described by Orla ndo et al.l (2006) in the context o f 



the study of supernova remnants (see also llVIiceli et al.ll2006l) . 
The results of numerical simulations are the evolution of tem- 
perature, density, and velocity of the plasma in one quadrant of 
the whole spatial domain. We reconstruct the 3-D spatial dis- 
tribution of these physical quantities in the whole spatial do- 
main, and we rotate the system about the z axis in such a way 
that the inclination of the binary orbit to the line-of-sight (LoS) 
is 35°, according to the values estimate d from radial velocity 
observations before the 2006 outburst ( Dobrzvcka & Kenvonl 



1994). Both the red giant and the white dwarf lie on the z axis, 
and such a configuration assumes that the 2006 explosion oc- 
cu rred at quadrature , in agree ment with the system ephemeris 
of lFekel et all d2000l) (see also lO'Brien et alj|2006l) . 

The emission measure in the j-th domain cell is emj = 



re^-Vj (where w^. is the hydrogen number density in the cell, 



Vj is the cell volume, and we assume fully ionized plasma). 
From the values of emission measure and temperature in the 
cell, we synthesize the corresponding X-ray spec trum, using 



the A strophysical Plasma Emission Code (APEC: ISmith et al 
200 lj) of hot collisiona lly ionized plasma. We a ssume so- 
lar metal abundances of iGrevesse & Sauvall (1998) (hereafter 
GS) for the CSM and abundances enhanced by xlO for the 
ejecta material (a s suggested by Chandra/HETG observations; 
Drak e et al . 2008). The spectral synthesis takes into account the 
thermal broadening of emission lines and the Doppler shift of 
lines due to the component of plasma velocity along the LoS. 

The X-ray spectrum from each cell is filtered through the 
absorption column density relative to the position of the cell 
in the spatial domain. The absorption consists of two com- 
ponents: 1) the photoelectric absorption by the ISM, assum- 
ing a column density A^h = 2.4 x 10 21 cm" 2 (according to 



of D op h = 1.6 kpc, e.g. iBode 1987 ) and 2) the local absorp- 
tion by the shocked CSM (with GS abundances) and by the 
ejecta (with GS abundances xlO) encountered within the blast 
wave. The ISM, CSM, and ejecta absorption components are 
computed using the absorption cross-sections as a function of 
wavelength from lBalucinska-Church & McCammon (1992). 

We integrate the absorbed X-ray spectra from the cells in 
the whole spatial domain. The resulting X-ray spectrum (orig- 
inating from the whole blast wave) is then folded through the 
Chandra/HETG instrument response, considering the source at 
a distance D op h = 1.6 kpc. The exposure time is assumed to be 
4xp = 10 ks (the same value of the Chandra/HETG observa 



tions of RS Oph analyzed by [Nelson et al.l 12008b iDrake et al 
2008). 



3. Results 

3.1. Hydrodynamic evolution 

In all the models examined, we found the typical evolution 
of radiative shocks propagating through an inhomogeneous 
medium: the fast expansion of the shock front with tempera- 
tures of several millions degrees and the development of dense 
and cold regions dominated by radiative cooling, as the for- 
ward and reverse shocks progress through the CSM and ejecta, 
respectively. As examples, Fig. [2] shows 2-D sections in the 
(x, z) plane of the distributions of mass density (on the left) 
and of temperature (on the right) for the models ND-E44-N7 
(i.e. without the EDE) and YD-E44-N7-L2 (with the EDE) at 
day 13.9. 

The radiative cooling is very effective in the inner portion 
of the remnant (dark regions in left panels of Fig. [2j made up 
mostly of ejecta material. Hot post-shock regions are domi- 
nated by thermal conduction that contrasts the radiative cool- 
ing and suppresses hydrodynamic instabilities t hat would de 



velop during the evolution of the blast wave (see lOrlando et al 



2005). In particular, the thermal conduction suppresses the 
Rayleigh-Taylor instabilities that would develop at the con- 
tact discontinuity and that would drive the mixing of ejecta 
material with the shocked CSM. In our simulations, the mass 
mixing mainly comes from ejecta evaporation driven by ther- 
mal conduction rather t han from hydrodynamic ablation (see 
Orlando et alj2005l[2008h . Fig.|2]also shows the distribution of 
ejecta within the blast wave: regions whose content is made up 
of original ejecta material by more than 90% (dashed contours 
in the figure) are mostly dense, cold and dominated by radia- 
tive cooling, whereas regions made up of a mixture of ejecta 
and CSM are hot and dominated by thermal conduction (see 
right panels of Fig. [2]). In all the cases, we found aspherical 
shock morphologies, with the blast wave propagating through 
inhomogeneous circumstellar gas distribution and being par- 
tially refracted around the red giant star. 

In models without EDE, the aspherical morphology is due 
to the fact that the origin of the blast wave is off-set from the 
origin of the wind density distribution by 1.5 AU. This off- 
set is very effective in collimating the blast wave and ejecta. 
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Figure 2. 2-D sections in the (x, z) plane of the mass density distribution, in log scale, (on the left) and of the temperature (on the 
right) in the simulations ND-E44-N7 (upper panels) and YD-E44-N7-L2 (lower panels) at day 13.9. The EDE lies in the (y,z) 
plane. The contours enclose zones whose content is made up of original ejecta material by more than 10% (dotted) and 90% 
(dashed). 



The shock front initially propagating away (toward) from the 
red giant goes through the r~ 2 (r 2 ) density gradient and, as a 
consequence, its propagation velocity and temperature are high 
(small): the blast wave appears elongated in the direction away 
from the red giant and the hottest plasma component emitting 
in the X-ray band lies to the side away from the red giant (see 
upper panels of Fig.[2j». The shock front propagating toward the 
red giant is strongly reflected by the high density medium and 
the resulting reflected shock contributes to sweep out the ejecta 
material away from the red giant (and from the white dwarf). 

In models including the EDE, the aspherical shock mor- 
phology is quite complex and originates from the propaga- 
tion of the shock through both the off-set red giant wind 
and the density enhancement at the equato rial plane (see also 
Girard & Willsonlll987tlWalder et alfeooil) . The latter compo- 



nent determines the collimation of the blast wave perpendicu- 
larly to the plane of the orbit of the central binary system and 
leads to a bipolar shock morphology distorted (by the off-set 
red giant wind) and converging on the side away from the red 
giant (see lower panels of Fig. [2]). In this case, the hottest X- 
ray emitting plasma component lies in the two poles to the side 
away from the red giant (see lower right panel of Fig. 12). A 
movie showing the evolution of 3D spatial distribution of mass 
density, in log scale, for the model YD-E44-N7-L2 is provided 
as on-line material: the point of view is along the y axis during 
the whole evolution, than it first rotates about the x axis and 
finally about the z axis up to an inclination of the binary orbit 
of 35° (i.e. the value estimated before the 2006 outburst). 
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Figure 3. Emission measure vs. tempera- 
ture distribution, EM(r), of the blast wave 
at day 13.9. Upper panels show the re- 
sult for models without EDE, lower pan- 
els for models with EDE. The black his- 
togram represents the average EM(T) esti- 
mated from t he analysis of Cha ndrafHETG 
observations jDrake et alj |2008). The factor 
10 44 5 assumes a distance for RS Oph of 1.6 
kpc. 



3.2. Emission measure vs. temperature 

From the model results, we derived the EM(T) distribution of 
the blast wave. As discussed in Sect. 12.21 from the distribution 
of mass density, we derive the emission measure in the j-th 
domain cell as emj = nj^Vj. The EM(T) distribution is then 
derived by binning the emission measure values into slots of 
temperature; the range of temperature [4 < log T(K) < 8] is 
divided into 15 bins, all equal on a logarithmic scale. Fig. [3] 
shows the EM(T) for models without EDE (upper panels), and 
with EDE (lower panels) at day 13.9, in the range of tempera- 
ture [5 < l og T(K) < 81, an d assuming a distance for RS Oph 
of 1 .6 kpc ( Bode et al]2007 ). The figure also shows the average 
EM(r) estimated from resonance lines o f abundant He- and H- 
like ions observed with Chandra/HETG dDrake et al.f 2008). 



In models without EDE, the shape of the EM(T) is char- 
acterized by a bump centered at temperatures between 1 and 5 
MK (see upper panels in Fig. [3). Increasing the initial energy of 
the explosion £bo from 10 43 to 10 44 erg or increasing the den- 
sity of the red giant wind at 1 AU from 2 X 10 7 to 2 X 10 10 cm -3 , 
increases the EM(T) distribution at all temperatures. However, 
in all the cases examined, the values of emission measure at 
temperatures around 10 MK are orders of magnitude lower than 



those derived from the observations (Dr ake et al.l 12008) and 



most of the shocked plasma is at temperatures below 5 MK. An 
initial energy of the explosion Ew > 10 45 erg (and an ejected 
mass M e j > 10~ 5 M Q ) and/or a density of the red giant wind at 
1 AU n w > 2 x 10 10 ctrT 3 (corresponding to a mass loss rate of 
> 5 x 10~ 6 M Q yr~' for a wind velocity of 30 km s -1 ) would 



probably reproduce the observed emission measure, leading to 
a bump around 10 MK in the EM(T) distribution. However, 
these values of Ebo and M e j are somew hat high for the tru e en- 
ergy and ejected mass of the explosion (1 Yaron et al . 2005) and 
the value of « w is also unrealistically large for the true density 
of the red giant wind. We conclude that models without the 
EDE fail to reproduce the values of emission measure derived 
from observations. 

In models including the EDE, the emission measure at tem- 
peratures above 5 MK is larger than that derived from the cor- 
responding models without EDE. The shocked plasma of the 
EDE leads to a bump around 10 MK in the EM(T) distribution 
(see lower panels in Fig. [3]). The emission measure of the bump 
depends on the initial energy of the explosion E^o- In models 
with E\,o = 10 43 erg, the EM around 10 MK is much lower than 
that derived from observations. Indeed the observed EM values 
can be reproduced with the models with E^o = 10 44 erg. The 
shape of the EM(T) distribution depends very slightly on the 
density of the red giant wind: decreasing n w , the EM slightly 
increases at temperatures above « 10 MK. On the other hand, 
the EM(T) distribution strongly depends on the thickness of the 
EDE: the model with the thinnest EDE predicts EM values at 
10 MK well below those observed. We found that model YD- 
E44-N7-L2 is our "best-fit" model for reproducing the average 
EM(T) estimated from observations (see lower right panel in 
Fig. El). 

Note that our models underestimate the emission measure 
values observed at temperatures > 20 MK. The "observed" 
EM(T) distribution at high temperatures is mainly derived from 
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Figure 4. X-ray light curves in the [0.6 - 12.4] keV band orig- 
inating from the nova outburst (solid line) and from ejecta 
heated by the reverse shock (dashed line) for YD-E44-N7-L2. 
The X-ray flux is obtained by integrating the synthetic X-ray 
emission in the whole spatial domain and takes into account 
the effects of local absorption by shocked CSM and ejecta. The 
reference dotted line is a power law with index a = -1.2. The 
X-ray luminosity is reported on the right axis assuming a dis- 
tance to RS Oph of 1.6 kpc. 



the analy sis of FeXXV (ytl .85), assuming GS abundances. 

d2008l) noted that EM values derived 



However, Drake 



from Fe ions appear to be anomalous and suggested that Fe 
probably suffers from effects of non-equilibrium ionization. In 
addition, as discussed in the next section, we found that a large 
contribution to the emission in the FeXXV (Al. 85) line comes 
from the shocked ejecta material, with GS a bundances xlO. As 
a consequence, the Fe XXV EM derived by Drake et al. (2008) 
would be overestimated. 

3.3. X-ray emission 

From the simulations, we synthesized the X-ray emission in 
the [0.6 - 12.4] keV band to match the observing range of 
Chandra/HETG, using the method outlined in Sect. 12.21 Fig. [4] 
shows the X-ray light curve in this band for YD-E44-N7-L2 
(i.e. our "best-fit" model), obtained by integrating the synthetic 
X-ray emission in the whole spatial domain and taking into 
account the effects of local absorption by shocked CSM and 
ejecta. The figure shows the X-ray flux, Fx, originating from 
the nova outburst together with the contribution to X-rays from 
ejecta heated by the reverse shock (dashed line). The X-ray flux 
reaches its maximum quite early, around day 4, and then de- 
cays as r 12 . The time of maximu m X-ray emission is consis- 
tent with Swift/XRT observations feode et al]|2006l) . although 
the observed X-ray flux is slightly higheid (by a factor « 2) 



than that derived from the simulation and decays as t~ (see 
also ISokoloski et al.ll2.006b instead of r L1 as derived from the 
simulation. The X-ray decay rate is mainly driven by radiative 
cooling, which has a significant impact on the post-shock tem- 
perature structure dO'Brien & Kahnll987l) . The contribution to 
X-ray luminosity from shocked ejecta is rather important dur- 
ing the early phases of the evolution (more than 80% for t < 5 
days) and then decreases down to 34% at day 15 due to signif- 
icant radiative cooling. 

Fig. [5] shows the map of X-ray emission integrated along 
the LoS and on pixels with size ~ 0.1 AU at day 13.9 for YD- 
E44-N7-L2. The map is drawn in linear (upper panel) and log- 
arithmic (lower panel) scales to highlight structures with very 
different emission levels. The plane of the orbit of the central 
binary system lies on the (y, z) plane and is assumed to be in- 
clined by 35° to the LoS, con sistent with the inclination o f the 
orbit from orbital solutions dDobrzvcka & Kenvonlll994t see 
also Sect. 12.2b . A movie showing the evolution of X-ray emis- 
sion for YD-E44-N7-L2 is also provided as on-line material: 
the dotted gray line in the movie marks the position of the for- 
ward shock. 

Fig.|6]shows a 3-D rendering of mass density, in log scale, 
at day 13.9 for YD-E44-N7-L2. We use two different color ta- 
bles to identify the plasma with temperature larger (red color 
table) or lower (blue) than 1 MK. The volume has been clipped 
to show the internal structure of the blast wave. The compari- 
son between Fig. [5] and Fig. fallows us to link X-ray emitting 
structures with plasma structures originating from the interac- 
tion of the blast wave with the inhomogeneous CSM. 

Most of the X-ray emission originates from an irregular jet- 
like structure, lying on the (y,z) plane, with a size of « 10 AU 
(see Fig . O ■ The structure is located to the side away from the 
red giant, and the X-ray emission is maximum in a region with 
size « 4 AU. By comparing the X-ray images in Fig. [5] with 
the 3-D rendering of mass density in Fig. [6] we note that the X- 
ray emitting structure corresponds to a region just behind the 
interaction front between the blast wave and the EDE (yellow 
region in Fig. [6}; the post-shock plasma at the equator is char- 
acterized by high density values (10~ 15 < p < 10~ 13 gm crrT 3 ) 
and temperatures ranging between 1 and 50 MK, with the high- 
est values of density and temperature to the side away from the 
red giant. 

The lower panel in Fig. [5] shows that faint X-ray emis- 
sion also arises from a double-lobed structure with a size of 
« 30 AU, surrounding the X-ray emitting jet. By compar- 
ing Fig. and Fig. [6] we found that the double-lobed struc- 
ture corresponds to regions downstream of the forward shock 
propagating perpendicularly to the plane of the binary orbit. 
These regions are characterized by density values ranging be- 
tween 10" 18 gm cm -3 and 10" 17 gm crrT 3 and temperatures be- 
tween 10 and 85 MK. If the orientation of the EDE is the same 



as discussed by lO'Brien et a l. (2006), we speculate that this 
double-lobed structure roughly corresponds to the synchrotron- 
emitting partial shell observed in high-resolution radio obser- 



However, note that the X-ray flux reported bv lBodeetall 12006) 
is corrected for the effects of absorption by intervening material, 



whereas the synthesized X-ray flux takes into account the absorption 
by shocked CSM and ejecta. 
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Figure 6. Three-dimensional rendering of mass density, in log scale, at day 13.9 for the model YD-E44-N7-L2. One quarter of 
the volume has been removed to show the internal structure of the blast wave. Red (blue) color table marks the plasma with 
temperature larger (lower) than 1 MK. The plane of the orbit of the central binary system lies on the (y, z) plane. 



vations dO'Brien et al.1 120061; iRupen et all 12008b although X- 
ray and radio brightness distributions may be not co-spatia{| 

3.3.1 . Line profile analysis 

We integrated the emission of the whole spatial domain and 
derived the spectra for the HEG and MEG traces as predicted to 
be observed with the Chandra/HETG, in order to compar e our 
anal ysis with the resu lts obtained by Nelson et a l. (2008) and 
by lDrake et alJ (120080 . As expected, the synthetic spectra show 
emission lines from different elements, forming over a wide 
range of plasma temperatures, and reflecting the broad nature 
of the plasma temperature distribution discussed in Sect. 13.21 



Following iDrake et al.l d2008l) . we analyzed the line pro- 
files of the most prominent spectral lines to investigate the 
origin of the broadening and asymmetri es revealed in th e 
Chandra/HET G observations (see , also. iNelson et alJ 120081) . 
Analogously to lDrake et alJ Q2008J), we restrict our analysis to 
HEG profiles, except in the case of O VIII (/118.97) which falls 
outside of the HEG range and is observed by the MEG. Table|2] 
summarizes the results of our analysis for the abundant He 
and H-like ions, together with the FeXVII (A 15.01) resonance 
line. Fig.|7]shows the line profiles for the abundant H-like ions 



3 Note that synchrotron radio brightness distributions also depend 
strongly on gra dients of ambient ma gnetic field strength or of ambient 



SiXIV 016.18), MgXII 018.42), and NeX 0112.13) observed 
by the HEG and O VIII (A 18. 97) observed by the MEG. Note 
that the line profiles are affected by the instrument profile that 
is known to be more effective for decreasing wavelengths. 

The line profiles exhibit broadening and asymmetries re- 
markably similar to those observed with Chandra/HETG. In 
velocity terms, most of the synthetic lines exhibit FWHM rang- 
ing between 1000 and 3000 km s and FWZI ranging be- 
tw een 2500 and 8000 km s~', in agreement with the results 
of INelson etaLI d2008l) and brake et all d2008b (and similar t o 
those observed in the IR bv lDas et al.ll2006l ; lEvans et al.ll2007h . 
The only exception is the bright He-like FeXXV (Al, 85) reso- 
nance line which exhibits much larger full-widths due to sub- 
stantial contribution to the emission from the other two compo- 
nents of the FeXXV triplet at A « 1.86 A and A » 1.87 Aand 
due to strong instrument effects. 

Drake et alJ 020081) noted that the observed line profiles are 



plasma density lOrlando et alj|2007t) . 



more peaked than expected for a spherically-symmetric shock 
and suggested that the emission is strongly enhanced in the di- 
rection of the plane of the sky. Our model predicts that most of 
the X-ray emission originates in a compact region propagating 
away from the red giant in the direction perpendicular to the 
LoS (see Fig. O. As a result, we found that the line profiles are 
peaked as those observed. 

In most of the lines, the centroids are systematically blue- 
shifted and the amount of the shift depends on the wavelength 
(see Fig.|7]and TablefJ): from the Si XIII (,16.64) line to the H- 
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Figure 5. X-ray images in the [0.6 - 12.4] keV band in lin- 
ear (upper panel) and logarithmic (lower panel) scales de- 
rived from model YD-E44-N7-L2 at day 13.9. The plane of 
the orbit of the central binary system lies on the (y,z) plane 
and is assumed to be incli ned by 35° to the line-of-sight 
( Dobrzvcka & Kenvonl 19941) . The dashed contour encloses the 
ejecta material leading to the largest photoelectric absorption. 



like O VIII (A 18. 97) doublet, the amount of blue-shift increases 
up to its maximum value. A similar trend has been found by 
Nelson et al.1 ([2008) from the analysis of Chandra/HETG and 
XMM-Newton/RGS spectra of RS Oph at day 13.9, although 
the amount of blue-shift found by these authors is systemati- 
cally higher than that predicted by our model. At variance with 
the other lines, the FeXXV (Al. 85) resonance line exhibits a 
significant red-shift. Again, this different behavior is due to 
contribution to the emission from the other two components 
of the FeXXV triplet. 




-4000 -2000 2000 4000 
Velocity [km s" 1 ] 

Figure 7. Synthetic velocity profiles for the H-like resonance 
lines of Si XIV, Mg XII, Ne X, and O VIII derived from model 
YD-E44-N7-L2 at day 13.9. Si, Mg, and Ne lines are observed 
by the HEG, O line by the MEG. The figure also shows the 
contribution of shocked ejecta material (red lines) and shocked 
CSM (blue lines) to the X-ray emission. 



As a result of the systematic blue-shift of emission lines, the 
line profiles tend to be more extended to the blue than the red 
(see Fig. Q and the line profile blue-shift and red-shift at zero 
intensity in Table |2]i with the only exception of the Fe XXV 
(Al. 85) line. As in the Chandra/HETG spectra, BSZI and RSZI 
are observed to decrease with increasing wavelength, and this 
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Table 2. Parameters (in units of [km s -1 ]) characterizing the 
profiles of the most prominent spectral lines in the HEG and 
MEG spectra" 





'ctr 


FWHVF 


FWZF 


BSZF 


RSZF 


FeXXV (21.85) 


1410 


11351 


29188 


13183 


16004 


S XVI (,14.72) 


-190 


3172 


8882 


4631 


4250 


SXV (25.03) 


-30 


2380 


6843 


3302 


3540 


Si XIV (26.18) 


-291 


2183 


5580 


2959 


2620 


Si XIII (26.64) 


-203 


1578 


4510 


2458 


2052 


Mg XII (28.42) 


-160 


1246 


3561 


1941 


1620 


Mg XI (29. 16) 


-196 


981 


2942 


1667 


1275 


NeX(212.13) 


-222 


988 


2471 


1458 


1013 


Fe XVII (215.01) 


-309 


499 


1297 


1008 


289 


OVIII (218.97) 


-355 


553 


1343 


1066 


276 



" Line profiles include the instrument profile that is known to be more 
effective for decreasing wavelengths; 

b Shift of the line centroid (negative values are for blue-shift); 

Full widths at half maximum; 

d Full widths at zero intensity; 

' Line profile blue-shift at zero intensity; 

f Line profile red-shift at zero intensity; 



is particularly evident in the O VIII (218.97) line profile in the 
MEG spectrum. 

The dashed contour in the lower panel in Fig.[5]encloses the 
ejecta material leading to the largest photoelectric absorptiorQ. 
Given the geometry of the system and the inclination of the or- 
bit, this material mostly absorbs the emission originating from 
the receding portion of the X-ray emitting jet-like structure (see 
also upper panel in Fig. ID that contributes to the red-shifted 
emission. In fact, the figure shows an area of attenuated X-ray 
emission located in the left lobe of the bipolar structure. The 
resulting absorption of this component explains the asymmetry 
and the systematic blue-shift of line profiles predicted by our 
model. 

We also used the tracer associated with the ejecta to deter- 
mine the contribution of shocked ejecta to the X-ray emission 
(see red lines in Fig. |7j. We found that the red-shifted emission 
of the ejecta component is, in general, more absorbed than that 
of shocked CSM: the asymmetry of observed line profiles is 
mainly due to the ejecta component. We also found that, in gen- 
eral, the smaller the wavelength the larger is the contribution of 
shocked ejecta to the emission; in particular, the shocked ejecta 
contribute « 60% to the emission of FeXXV (21.85) line and 
a* 40% to the emission of Fe XVII (215.01) line. The different 
contributions of ejecta to X-ray emission (together with pos- 
sible deviations from equilibrium of ionization) coul d explain 



the an omalous EM values derived from Fe ions by lDrake et al 



d2008l) . Among the lines analyzed, that with the largest contri- 
bution of shocked ejecta (more than 70%) is the S XV (25.03) 
line. 

4 Eq. fallows us to trace the ejecta material in the CSM. 



4. Summary and conclusions 

We have investigated through detailed hydrodynamic modeling 
the origin of asymmetries and broadening of the emission lines 
observed with Chandra/HETG during the 2006 outburst of RS 
Oph. To our knowledge, these simulations represent the first 
attempt to model the 3D structure of the blast wave originat- 
ing from a nova outburst that simultaneously considers radia- 
tive cooling and thermal conduction (including heat flux satu- 
ration). Our findings have significant implications on the diag- 
nostics of the blast wave during the early phases of evolution 
and lead to several useful conclusions: 



1. 



In all the cases, we found that the nova remnant is highly 
aspherical, with the morphology of the blast wave deeply 
influenced by the pre-existing inhomogeneous CSM. Even 
in models without the EDE, the off-set of the nova explo- 
sion from the origin of the wind density distribution causes 
an effective collimation of the blast wave and ejecta. 
Models without the EDE cannot reproduce the values of 
emission measure derived from Chandra/HETG observa- 
tions unless the outburst energy is £t,o ^ 10 45 erg (and the 
ejected mass in the explosion M e j » 10~ 5 M ) and/or the 
density of the red giant wind is « w » 10 10 cuT 3 at 1 AU; 
however, for the physical parameters characterizing RS 
Oph, these values seem too high for the true outburst en ergy 
and ejected mass in the explosion dYaron et al. 2005h and 
for th e density of the red giant wind dDobrzvcka & Kenvon 



1994). Models including the EDE reproduce in a natural 
way, without any further ad hoc assumption, the values of 
emission measure derived from Chandra/HETG observa- 
tions if E w ~ 10 44 erg and M ej * 10~ 6 M Q . 
Our "best-fit" model (YD-E44-N7-L2) predicts that, at day 
13.9, most of the X-ray emission originates in a region with 
size « 4 AU localized at the interaction front between the 
blast wave and the EDE to the side away from the red gi- 
ant; such an X-ray emitting region propagates in the direc- 
tion perpendicular to the LoS. As a result, the synthetic line 
profiles are more peaked than expected for a spherically- 
symmetric shock in nice agreement with the observations. 
The synthetic line profiles are asymmetric and slightly 
blue-shifted and they are remarkably similar to those ob- 
served. We found that the observed asymmetries are due 
to substantial X-ray absorption of red-shifted emission by 



ejecta material, confirming the conclusion of Drak e et al 
(2008). 



5. Both shocked CSM and shocked ejecta contribute to the 
observed X-ray emission. The asymmetry and blue-shift of 
emission lines are mainly due to the shocked ejecta com- 
ponent which is more affected by X-ray absorption of red- 
shifted emission. In general, the contribution of shocked 
ejecta to the X-ray emission decreases with increasing 
wavelength. 

Our model shows that the broad range in the plasma tem- 
perature and the asymmetries of line profiles observed in 
Chandra/HETG spectra are due to the interaction of the blast 
wave with the pre-existing inhomogeneous CSM. In addition, 
the asymmetric nature of the CSM into which the early blast 



S. Orlando et al.: 3D modeling of the asymmetric blast wave from RS Ophiuchi 



11 



wave is driven is also responsible of the apparent shock colli- 
mation in the plane of the sky. Note that the jet-like ejection 
emitting in the X-ray band predicted by our model does not 
coincide with the synchrotron jet obse rved in the radio band 
dO'Brien et al.l2006l;lRupen et al.ll2008l) : the former propagates 
in the plane of the orbit of the central binary system, whereas 
the latter propagates perpendicularly to the same plane. 

Our analysis confirms that the comparison of high quality 
data with detailed hydrodynamic modeling can be a powerful 
tool to study and diagnose the physical properties of outbursts 
from recurrent novae and to interpret the observations. In par- 
ticular, our results may provide useful constraints on the cir- 
cumstellar gas distribution of RS Oph and suggest the existence 
of a dense region of the red giant wind at the equatorial plane. 
We also show that most of the early X-ray emission arises from 
this region as the blast wave from the 2006 outburst deceler- 
ates. 

The scenario of an equatorial density enhancement is con- 
sistent with observations of the 2006 outburst in other wave- 
length bands, for instance VLBA radio synchrotron observa- 
tions of an evolving ring-like structure (O'Brie n et al. 2006; 
Rup en et aL 20081) and HS T observations of a bipolar nebular 
structure dBode et al.ll2007 ). Also, the existence of the equato- 
rial density enhancement may explain the presence of dust in 
the circumstellar environment in wh ich the explosion occurs 
dEvans et al. 1120071: iBarrv et alj|2008l) . 

Our findings are also consistent with the asymmetric blast 
wave evolution described analytically by iGirard & Willsonl 
(119871) and with the results of hydrodynamic models describ- 
ing the evolution of the circumstellar environment in b inary 
systems (Mas trodemos & Morris! 1999llWalder et al.ll2008l) . In 
particular, in systems comprising a red giant star, these models 
predict an equatorial density enhancement created by a spiral 
shock wave caused by the motion of the stars through the cool 
wind of the red giant; the enhancement is most pronounced in 
systems with the smallest binary separation, comparable to that 
in RS Oph. 



It is worth noting that, in Wald er et alJ (12008), the blast 



wave appears more collimated perpendicularly to the plane of 
the orbit of the central binary system than in our simulations. 
The reason is that the CSM immediately around the white 
dwarf (within 1/3 AU from the dwarf) is characterized by a 
rather co mplex density struc ture (see the lower right panel of 
Fig. 2 of Walder et al. 2008) determining the strong collima- 
tion of the blast wave, at the very beginning of the nova erup- 
tion (not described by our model). Since we assume the initial 
blast wave to be spherical (see Sect. 12. IK our model does not 
take into account such a strong initial collimation of the shock. 
To evaluate the effect of the initial condition on the evolution 
of the blast wave, we compared our "best-fit" model with an 
additional simulation with an identical setup but assuming an 
initial ellipsoid- shaped blast with the major axis perpendicu- 
lar to the orbital plane. Such an initial condition describes the 
strong shock collimation occurred during the early phases of 
the nova eruption. Taking into account this initial shock colli- 
mation, we found that the shock expands very rapidly toward 
the perpendicular to the orbital plane and, after few days ofeyo- 
lution, its shape is similar to that found bv lWalder et alJ (12008b . 



Nevertheless, we also found that the main features characteriz- 
ing the X-ray emission (broadening, asymmetries, blue-shift of 
emission lines) arising from the blast wave do not change qual- 
itatively. Therefore, the main conclusions of this paper on the 
origin of the X-ray emission and of asymmetries and broaden- 
ing of emission lines observed with ChandrafHETG remain. 

The comparison of our model results with Chandra/HETG 
data suggests that the mass of ejecta in the 2006 outburst was of 
the order of 10~ 6 M . From the analysis of the spectra collected 



with th e Rossi X-ray Timing Explorer (RXTE). ISokoloski et al 



(2006) derived the ejecta mass to be of a few times 10~ 7 M B on 
the basis of some assumptions (time of transition to the Sedov- 
Taylor phase, shock speed during the ejecta-dominated phase, 
mass density inside the binary). However, taking into account 
the uncertainties in the determination of the se values, the mass 
of ejecta derived bv lSokoloski et al.l d2006l) can range between 
few times 10~ 7 M and few times 10~ 6 M , i.e. a range of val- 
ues consistent with that derived by our model. 

Considering the mass of ejecta in the 2006 outburst of the 
order of 10~ 6 M , it comes out that the white dwarf is increas- 
ing in mass if its growth rate is larger than 5 x 10~ 8 M yr _1 . 
For lower values of the growth rate, the nova would throw off 
more mass than the white dwarf may have accre ted in the inter- 
vening 22 years. Recently lHachisu et alJ (120071) have estimated 
the white dwarf mass in RS Oph to be 1.35 + 0.01 M and 
its growth rate t o be about 10~ 7 M yr~' in average (see also 
Kato et al. I l2008l) . According to our results, therefore, we con- 
clude that the white dwarf mass is effectively growing up and 
RS Oph could be the progenitor of a SN la as the white dwarf 
reaches the Chandrasekhar limit. 

It will be interesting to expand the present study, including 
the ambient magnetic field, to investigate the evolution of the 
blast wave in later evolutionary phases and to make predictions 
on the synchrotron radio emission. The detailed comparison of 
model results with observations may lead to a major advance in 
the study of interactions between the blast wave and the mag- 
netized CSM in recurrent novae, and may provide important 
insight into SNe la, of which recurrent novae are believed to be 
progenitors. 
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